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ABSTRACT 


This paper describes the design of the robust PI/PID controller for the higher 
order interval system via its reduced order model using the differential 
evolution (DE) algorithm. A stable reduced interval model is generated from 
a higher order interval system using the DE in order to minimize the cost and 
reduce the complexity of the system. This reduced order interval numerator 
and denominator polynomials are determined by minimizing the Integral 
Squared Error (ISE) using the DE. Then, using reduced order interval model, 
a robust PI/PID controller is designed based on the stability conditions for 
determining robust stability of interval system. Finally, using these stability 
conditions, a set of inequalities in terms of controller parameters is obtained 
from the reduced order closed loop characteristic polynomial. Then these 
inequalities are solved to obtain robust controller parameters with the help of 


a DE algorithm. The designed, robust controller from the reduced order 
interval model will be attributed to the higher order interval system. 
The designed PI/PID controller from our proposed method not only 
stabilizes the reduced order model, but also stabilizes the original higher 
order system. The viability of the proposed methodology is illustrated 
through the numerical example of its successful implementation. 
The efficacy of the proposed methodology is also evaluated against 
the available approaches presented in the literature and the results were 
successfully implemented. 


Robust controller 
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1. INTRODUCTION 

In most of control engineering problems, design and development of a controller are an important 
part among the researchers, because the simulation and design of a controller of the higher order system is 
not an easy task. The cost and complexity of the controller increase with system order. To conquer this 
problem of simulation, analysis and design of a controller for the higher order system, a reduced order 
model is essential for the original higher order system. In the closed loop system, a controller designed for 
this low order reduced model, will also stabilize the original higher order system. Hence, from the last three 
decades, model order reduction of the higher order system is an important aspect for researchers. For this 
reason, many methods are available for lower order models which are expected to approximate 
the performance of a higher order of linear time invariant system. Various model order techniques have 
been proposed for the order reduction of higher order continuous time as well as discrete-time systems. 
Among these methods, the familiar and important methods are discussed here. An aggregation method was 
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proposed control the large scale dynamic system [1]. Pade approximate method [2] has very useful features 
for research due to being computationally simple. But the approximant obtained by this method often leads 
to be unstable even though the original system is stable. In view to conquer the problem of stability, 
the Routh approximation method [3] was introduced. But, the Routh approximation method fails to give 
well approximate at the low frequency level of the higher-order system. Several methods [4]-[7] have been 
given to improve the Routh approximation method for obtaining a reduced order model. All these methods 
are available for systems with fixed coefficients only. However, many control engineering problems like 
flexible manipulator system or nuclear reactor system represented by their mathematical model in a wide 
range of unknown operating conditions are bound under a certain amount of uncertainty. These uncertain 
parameters can be represented by bounded interval. This motivates researchers for an invention of 
the classical techniques for interval systems, or even for more general uncertain systems. Since last decade, 
much effort has been made in the design, simulation and control of plants with parametric uncertainty. 
Unfortunately, there are few methods available in the literature for the design of plants with parametric 
uncertainty through model order reduction in the entire range of operation. Some of the model reduction 
methods for fixed co-efficients mentioned above are extended to order the reduction of higher order 
interval systems through interval arithmetic to derive lower order approximates. Routh approximation 
technique [8] based on the direct truncation of the Routh table for interval systems using interval arithmetic 
and extending the Routh-Pade approximation reduction technique [4] to the interval system. Hwang 
et al. [9] pointed out that this method has sometimes failed to generate stable reduced order interval models 
for stable original systems. To improve the effectiveness of order reduction of interval systems, many 
mixed methods are proposed in [10]-[12]. The stability analysis of interval polynomial was presented in 
[13] using Hurwitz polynomial. 

In the area of control engineering, a very less work has been made in the design of controller for 
higher order systems using model order reduction technique. Hence, it is essential for researchers to design 
a controller for the higher order system through the reduced order model. Here few methods are introduced 
in the design of controllers which are available in the literature. From Lower order generalized aggregated 
model [14], a stabilizing controller is designed for the reduced order model, if applied to its original higher 
order system; it does not always guarantee the stability of its closed loop system. A lower order control 
technique proposed [15] using a model order reduction of an optimal closed loop system. But from this 
control technique, the higher order system stability not retained for a stable original system. Most of 
research work has carried out in this direction and a few methods are available. A controller scheme [16] is 
proposed, which is restricted to particular Davison reduced order model only and requires the system states 
to be available for feedback. Certain stability conditions are mentioned in [17] that the controller designed 
is not only stabilizes the reduced order model but also stabilizes the original higher order system. He also 
showed that the approximate inclusion of any unstable real modes of the higher order system in the lower 
order model will retain the existence of such simultaneous stabilizing controller. A special care for any 
model reduction scheme [18] has taken to obtain stabilizing controllers from a reduced order model. 
In the direction of Arno Linnemann, a dominant pole retention technique [19] is presented in a controller 
design for continuous time systems through reduced order model. All these available methods are used only 
for fixed coefficients of a polynomial/system. 

In control engineering applications, tuning of PI, PID and lag/lead controllers have been widely 
used in industries for several decades. However, many important results have been recently reported on 
computation of all stabilizing P, PI and PID controllers in [20]-[22]. These controllers are widely used in 
various process control applications. The controller performance was compared based on settling time, 
percent of overshoot and stability analysis of a given system with fixed coefficients. But, many control 
engineering problems represented by their mathematical model in a wide range of unknown operating 
conditions are bound under a certain amount of uncertainty. The large uncertainty present in the control 
system causes degradation of system performance and destabilization. The model is known approximately 
and hence it is necessary to incorporate the robust in design. Therefore, robust control present in these 
uncertainties is very important for plant operation under stabilized condition. This necessitated a robust 
controller design which could stabilize the plant for all the operating conditions. Hence designing a robust 
controller for parametric uncertain plants having unknown, but bounded parameter uncertainties has 
become the problem of research nowadays. To minimize the stated uncertainties, many solutions were 
proposed in the literature for the simulation, design and tuning of controllers [23]-[24]. Recently, affordable 
results have been reported on computation of all stabilizing P, PI and PID controllers. Therefore, after the 
renowned theorem of Kharitonov [13], the stability analysis of polynomials due to parameter uncertainty 
has an important aspects of researchers. As per Kharitonov, the interval polynomial which assesses robust 
stability conditions that four specially constructed extreme polynomials are Hurwitz. Robust stability of 
interval polynomial is also discussed by many researchers’. In [25], necessary and sufficient stability 
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conditions for interval polynomials are proposed using the results of [26] for fixed polynomials. 

In this paper, model order reduction of interval systems is carried out by using the Kharitonov’s 
polynomial [13] and differential evolution [27] using the ISE method. The numerator and denominator 
polynomials of the reduced model are obtained by minimizing integral squared error between the transient 
response of the original higher order system and the reduced order model pertaining to a unit step input. 
Thus, the stability is guaranteed for the reduced order system if the original higher order system is stable 
and the responses matching between original higher order system and the reduced order model. Then 
a PI/PID controller is designed for reduced lower order interval process plant based on the necessary and 
sufficient stability conditions [28]. These conditions are used to derive a set of inequalities in terms of 
controller parameters. The inequality constraints in the polynomial are solved consequently to obtain 
the controller parameters with the help of DE. The PI/PID controllers designed in this proposed method not 
only stabilizes the reduced order model, but also the original higher order system. The efficacy of 
the proposed method is demonstrated by implementing with a typical numerical example available in 
the literature. In comparison with the methods available in the literature [23], the proposed method in this 
paper is simple and involves less computational complexity. This paper is organized as follows: Section 2 
gives robust of stability conditions of interval polynomial. Section 3 presents the problem formulation and 
the proposed order reduction algorithm. Section 4 proposes the design procedure for robust stabilizing 
PI/PID controller in section 5 proposed method is applied to design a robust PI/PID controller for higher 
order interval process plant with a suitable numerical example and the final conclusion is given in section 6. 


2. ROBUST STABILITY CONDITIONS OF INTERVAL POLYNOMIAL 
Consider an interval polynomial ‘P(s)’ of order n of the form given below. 


P(s) = pps” + Pris! Epa sa ps! +... + PIS + Po, 


Where Pi £! 4-7! for i=0,1,2,3,.....,n. 

According to Anderson, Jury and Mansour [29] the necessary and sufficient condition for robust 
stability of interval polynomials of order n<3 is positive lower bounds on the coefficients of an interval 
polynomial. Therefore, consider an interval polynomial of order n=1 


P(s) = Xl opis, where p; €[a;,b;]. 
P(s) = pS + Po = [a;,b; Is + [ag. bo]. 


Therefore, as per Anderson et al., [29], the robust stability condition is 


aj >Oandag >0 1.e. a;>0 for 1=0,1. 
Similarly for order n=2 
P(s) = X- opis = pos? +pıSs + po = [az, b2 ]s? +[a4, b1 |s+[ag, bo]. 
Therefore, the robust stability condition is 
ay >0, ay >0 and aọ>0 Le. a; >0 for 1=0,1,2. 
Lemma 2.1 Consider a real Hurwitz polynomial Q(s) of the form (1) 
Q(s)= qns” kir Aage +....+q1Ss+q0 (1) 
is OF eee ee n 


Where q; is real and positive, qq >0. 


If any complex number Z such that R, >0, f (z)| > f (—z)|, moreover, 





f(z) > |f(—2)| , where C is a Closed 
zon C zon C 


contour, then, according to Routhe’s theorem [28] the following two polynomials can be formulated (2-3). 
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1 
Qo = 51 Qs) +A- (2) 


X 


Qi = Q9- 2 (3) 
2s 


X 


Theorem 2.1:For stability of Qs)the two polynomials Qg and Q,formed by the alternate 


coefficients of a Hurwitz polynomial in accordance with equations (2) and (3) must have negative real 
zeros. The proof of this is given in [28]. 


2.1. Necessary conditions for the stability of an interval polynomials 
Consider an interval polynomial of order n > 3 of the form (4) 


P(s) = Pps” tay gg +...+ pisi +...+ P1S + PQ; (4) 


Where p; €[a;,b;] for i=0,1,2,3,......n. 
The necessary conditions for an interval polynomial to be stable is given as (5) 


b; 2a; >0 fori= 0,1,2,3,...... n. (5) 


2.2. Sufficient conditions for stability of interval polynomial for n > 3 
A real coefficient interval polynomial of degree n is of the form (6) 


P(s) = pns? +pp_iS” | +.....+ PIS + PO (6) 
Where n> 4, 

Po € 140,01, pı € [47,67], p2 € [a2,b2 ],p3 € [43,03 J... Pp EL Ay rby I 

Theorem 2.2: A real coefficient interval polynomial of degree" 24 given in equation (6) is said to 
be robustly stable if and only if the following conditions are satisfied. Therefore, the necessary and 


sufficient conditions for stability of interval polynomials of order n>3 are derived from [28] and they are 
presented in Table 1. 


Table 1. Robust Stability Conditions for Various Higher Order Interval Polynomials 


Order of thePolynomial Necessary Conditions Sufficient conditions 
n 
3 a a? > 3bọb2 
Where i = 0,1,2,3. 
S a 
4 2 di >0 aî > 4bob4 and Ego 
0 b3 
Where i = 0, 1,2,3,4. 
n 
5 aon a2 > 4bob4 and a3 > 4bjbs 
Where i = 0, 1,..4, 5. 
n 
6 žaj >0 
0 a3 > 3bọb4 and a4 > 4bjbs. 
Where i = 0, 1,., 5,6. 
5 0 
7 2a; á a3 >3bob4 and aĝ > 3bybs. 
Where i = 0, 1,., 6, . 
< a 
a; >0 2 2 —a?2 
8 2 i a >3bıb5, b4 > 4agag and 5 = <0. 
Where i = 0, 1,,.., 7,8. 
n 
: —a —a 
9 ee b > 4agag bs > 4ajag,—* <0 and aoe 
6 7 
Where i = 01,,..,8,9. 
n 
; —a 
10 2a; >0 bi > 4agag ‘ be > 4ajag, bg > 4a7a 9 and = <0 
7 


Where i =, 1, 9, 10. 
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3. Proposed order reduction algorithm 
3.1. Problem formulation 
Consider a higher order continuous time interval system given by the transfer function (7): 


= N(s,C) [C5 56; ]1+1C7,C} Is Pitti Cy Cm is 


G(s,C,D) = — + — + = + 
D(s,D) [Do ,Do 1+[D] D] ]1s+.....+[ Dn, Dn Is" m<n 





(7) 


Ce = N) C 
Where [Ci >Si Ifor 1=012..--m are numerator coefficients of SGD) withG and “i as lower and upper 
+ 
] 


: C; ,C:; „D? i 
bounds of interval [Sis Ci i Viie r=0, L2 


; D; ‘ AE 
respectively, and ea ™ are denominator coefficients of 


+ — p+ 

and A as lower and upper bounds of interval Cii respectively. 

It is desired to synthesize a robust controller for this higher order interval plant. This can be 
achieved by approximate the higher order plant into the reduced order model using the proposed model 
order reduction algorithm and a PI/PID robust controller is designed for the reduced order model using 
the proposed algorithm given in section 4. Finally the designed controller is attributed to the higher order 
system for robust stability. 

The k" order reduced interval model is obtained from our proposed method is given as (8) 


G(s, C, D) with D; 


z = 3 k-1 
N(s,c) [co cg] + [cy c] Is + ..-----+ [ek-1 Ck ]8 


G,. (s,c,d) = = =e =a ane 
D(s, d) [do, do] + [dy , d} |S +...... + [dy, dy |s 





(8) 


— + = + 
ae: "i are alee c; 
Where [°° / for i=0,12...-k—1 are numerator coefficients of (5:°4) withĉi and “i as lower and upper 


7 q? ; . 8 
] respectively, and [di di I for i=0.12..-k are denominator coefficients of 
= k , Z JF 
G(s,cd) with “i and “I as lower and upper bounds of interval bay eed 


er cu 
bounds of interval EG 


respectively. 


3.2. Order Reduction Of Interval System 
According to the Kharitonov’s theorem, the interval system can be represented into four fixed 
parameter Kharitonov transfer functions. They are given as (9-12): 





Gls. cp ya MO) Go +O s+ Es" + ni Lagat Case 

D'(s,D) D) +D,s eae ee Dis (9) 
Puen Noe) Co tCrst Gs" + Paes: Cas 

D*(s,D) D} +Dj's+D3s? +......D, s” (10) 
isop- _ So tC s+ as" + ee Chs” 

D°(s,D) DË + Dy s+ Dys? +......Dp s" (1) 
Fisch N4(3,C) _ Cy E +078" aa 

D*(s,D) Dj + Dj's + Dzs’ + ae Sa eve D; s” (12) 

The above four Kharitonov’s transfer functions are, in general represented as (13 
g p 
i 

Gi(s,C,p)2=2 CA 

D (sD) for i=1,2,3,4 (13) 


These higher order four fixed parameter Kharitonov’s transfer functions are approximated into 
the reduced order model using following procedure. Let the four fixed parameter reduced order models of 
the above equations (9-13) obtained by the proposed method are defined as (14-17) 
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T +c, 8+c48" — k-l 
G} (s,c E ea a E 
9s 1 m = ; = - 
d,(s,d) dg +d, s+d3s +...... d,s a 
mie Sets - k- 
ep a AG ES 
9 9 2 -— 5 = 7 
d<(s,d) do + dj's + d3s tod S 
CEE SFOS + k- 
G% (s.c jy re) 6p T ee 
3 9 3 = = 5 7 
d? (s,d) do +d s+d7s^ +... dgs T 
4 = n4 (s,c) = on +¢,'s +0987 INVERSE cis 
OOR ea a a ee E 
d; (s,d) dg tdi sdas +...... dgs T 


The above Four Kharitonov’s transfer functions are, in general represented as (18): 


. ni S,C 
G} (s,c,d) = kl ) 
d! (s 

k ` 





°) For i=l. 2.3.4. (18) 


In this work, the numerator and denominator coefficients of the reduced model are obtained by 
minimizing the objective function ’J‘, which is the error between the original higher order system and 
the reduced order system. Therefore, it is represented in the form (19): 


_ foo _ Jk 2 
J=j YW- y dt (19) 


Mathematically the Integral Square error can be represented as (20) 


M 
J=5lyt)- y“ Wf 
0 (20) 


Where, y (t) is the unit step response of higher order and y* (t) is the unit step response lower order 
system at the instant in the time interval 0 < t < M, where M is to be chosen. The objective is to obtain 
a reduced order model, which is closely approximate original system. The objective function is to minimize 
ISE by using DE. 


3.3. Proposed DE algorithm for order reduction 

DE 1s a stochastic, population based direct search optimization algorithm introduced by Storn and 
Price in 1997 [27]. DE works with two populations; old generation and new generation of the same 
population. NP is the size of the population and it is adjusted. The population consists of real valued vectors 
with a dimension D that equals the number of design parameters/control variables. The population is 
randomly initialized within the initial parameter bounds. The three main operations carry optimization 
processes are: mutation, crossover and selection. In each generation, individuals of the current population 
become target vectors. For each target vector, the mutation operation produces a mutant vector. 
The crossover operation generates a new vector, called trial vectors, by mixing the parameters of the mutant 
vector with those of the target vector. If the trial vector obtains a better fitness value than the target vector, 
then the trial vector replaces the target vector in the next generation. 


3.3.1. Initialization 
Define upper and lower bounds for each parameter of the reduced order model (21) 


L U 
Xj SA Ajy (21) 


Eoy 
A 


Randomly select the initial parameter values uniformly on the intervals ‘and the elements of each 


xj 
individual +° are given by (22) 
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j _ yj j j 


Where 1 = 1,2,...,N and j= 1, 2,...,D. 
Where N is the population size, rand (0,1) is a random number uniformly distributed between 0 and 1, D is 
the number of control variables. 


3.3.2. Mutation 

Mutation expands the search space. DE undergoes mutation operation after initialization. 
In mutation operation, it produces mutant vector Vi,c, with respective to each individual Xi,c, so called 
target vector, in the current population via mutation strategy (23): 


VG =Xi,G t+ Fl Xpest.G —Xi,g J+ F(Xr1G — X712,G )- (23) 


i X; X X 
For a given parameter vector “iC two vectors "tO and “>C are selected randomly such that 


the indices r1, r2 is distinct. The mutation factor F is a constant from [0, 2] is called the donor vector. 


3.3.3. Crossover 
Crossover incorporates successful solutions from the previous generation. After mutation, DE 


. U; . X: G, 
undergoes crossover. The trial vector >S is developed from the elements of the target vector, +S? and 
Vi G: 
the elements of the donor vector, +S` (24) 
i = yar rand ;(0,1)s CR Jor j= Jrand )» 
= x! otherwise 
»§ (24) 


Elements of the donor vector enter the trial vector with probability CR (cross over rate) set to [0; 1]. 


3.3.4. Selection 
The newly generated values of trail vectors exceed the corresponding upper and lower bounds; we 
initialize them randomly and uniformly within the pre-specified range (25): 


j -Ul ; j 
Xi G =U; gif fU; G) <f(Xi a) 


_ wyj . 
= Xi g otherwise 


(25) 


The trail vector Xis is compared with trail vector Uis and the one with lowest function value is 
admitted to the next generation. Therefore the four k" order reduced Kharitonov’s transfer function 
denominators and the numerators are obtained by minimizing integral square error using Differential 
Evolution Algorithm. 

Finally the reduced order interval model is obtained by the following equation (26): 


k-l f 

2 [min 6j;),max jj ls? 
GÈ (s,c,d) = E i 
ps [min@;;),max jj) Js? 


icv For i=1, 2, 3,4. (26) 


For this reduced order model, a PI/PID robust controller is designed based on the minimization of 


J =(Kp - Kp)” + (Ky -Ky)* q J=(Kp-Kp)* + (Ky -Kp)* +(Kp - Kp)” 


the objective function for PI controller an 


for PID controller using the same DE algorithm. 
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4. Design Procedure for Robust stabilizing PI/PID Controller 


-Ge(S) 


Let the stabilizing controlle is considered to be a PI/PID controller transfer function of 


the form given below (27-28). 





G.(s)=Kp PAL = Se 

i c(S) for PI controller (27) 
G,(s) = Kp + aE + Kps = ae 

à c(8) for PID controller (28) 


Where ŚP = Proportional gain, “! =Integral gain and *P derivative gain. 
Now the system with robust stabilizing controller for parametric uncertainty is as shown in Figure 1. 


Robust Controller 





Figure 1. Block diagram of closed loop system with PI/PID controller 


Then the closed loop transfer function with a PI / PID controller can be defined as (29) 


Tís)= G.(s)G,(s,c,d) _ N-(5)N,(s,c) 


~1+G,(5)G,(s,6,d) Ne(s)N;(s,c)+De(s)D,(s,d) (29) 


Where Nr and Prs-® are the numerator and denominator polynomials of the reduced order 


Ne(s) De (s) G,(s) 


interval plant, and are the numerator and denominator polynomials of respectively. This 


PI/PID controller robustly stabilizes the interval plants family, if for all ceC and¢d<D then 
the characteristic polynomial of a closed loop transfer function given in equation (29) has all zeros have 
negative real values. Now apply the necessary and sufficient conditions of robust stability conditions 


given in Table.l to the closed-loop characteristic polynomial NeSNr(.9+De(SD59 Which 
leads to a set of constraints in terms of controller parameters. Then these constraints are solved 


jo(Kp=Kpy + (K1 - K9)? 


by using DE so as to minimize the objective function for PI controller and 


_ 0,2 0,2 0.2 
J=(Kp—Kp)" +(K1-KI) +(Kp-KD) for PID controller to obtain controller parameters. Then after 


obtaining the controller parameters, form four sets of Kharitonov’s polynomials [13] to check the stability 
and the closed-loop step response to verify the results. 

Application of the DE algorithm for determining the controller parameters is as follows: 
Step 1: Initialization: 

Define upper and lower bounds for each controller parameter of the PI/PID controller (30) 


L U 
Xj <XjilSXj. (30) 


The initial parameter values are randomly selected uniformly within the allowable range (31) 


LSU j _ yj j j 
[X; Xj ], Xio SXi + rand (X mx =A (31) 
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Where I = 1,2,...,N and j= 1, 2, DN is the population size, rand (0; 1) is a random number 
uniformly distributed between O and 1, D is the number of controller parameters. The minimum and 


maximum values of the control variables for PI/PID Controller i.e. $P ; Ki and ŚP are selected as 
0< Kp <10 and < KI <5 for PI controller. 


0< Kp <10, 0< Ky <5 ang 9<Kp <5 forPID controller. 


Step 2: Mutation 

DE undergoes mutation operation after initialization. In mutation operation, it produces mutant 
vector Vic, with respective to each individual Xi, so called target vector, in the current population via 
mutation strategy (32): 


Vig =Xig tl  Xpest.Gg ~ Xig )+ l(X+1gG — X124 >) (32) 


For a given parameter vector “+S two vectors “HS and ">S are selected randomly such that theindices 


rl, r2 is distinct. The mutation factor F=1.2 is chosen from [0; 2] is called the donor vector. 
Step 3: Crossover 


U; , 
bG is developed from the elements of 


G 3 


After mutation, DE undergoes crossover. The trial vector 


Xi G 


V: 
the target vector, and the elements of the donor vector (33), ” 


ulg = v? Gif ( rand ;(0,1)< CR Jor j = jrand ) 


= ye otherwise 


(33) 
Elements of the donor vector enter the trial vector with probability CR = 0.6 set from [0; 1]. 


Step 4: Selection 
The newly generated values of trail vectors exceed the corresponding upper and lower bounds; 
we initialize them randomly and uniformly within the pre-specified range (34): 


XJ Gq =U/ pif FUjfGgIES F(X Gg) 
= x) otherwise 


g (34) 


The trail vector Xi is compared with trail vector Uig and the one with lowest function value is 
admitted to the next generation. Finally the results of the best parameters of a PI/PID controller and fitness 
are obtained 


5. NUMERICAL EXAMPLE 
Consider a wing aircraft [23] whose transfer function with parametric uncertainty is given (35) 


[54,74]s+[90,166] 
St ar SE Ia IRE Se See Dae Ses 
S +[2.8,4.6 ]s~ + [504,808 ]s~ +[30.1,33.9]s+[-0.1,0.1] (35) 


The closed loop transfer function of the above interval system is given by (36) 


[54 74]s+[90 166] 
[1 1]sf +[2.8 4.6]s° +[504 808]s* +[841 1079]s+ 
[89.9 1661] (36) 


T(s)= 


This higher order closed loop transfer function of interval system can be represented as four fixed 
parameter Kharitonov transfer functions that are given as (37-40): 
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54s +90 
L ERTE a ETTE 
Ss +4.6s +808s° + 84.1s +89.9 


74s +90 

E ae oe ETOR ET 

Ss +2.85 +80.8s° +1079s +89.9 
54s +166 

ae E a a aE 

s` +4.6s° +504s° +&64154+ 1661 
745+ 166 

iS a are ee eee 

s +2.85 +50.4s° +1079s + 1661 


O 241 


(37) 


(38) 


(39) 


(40) 


It is difficult to analyze the higher order system, because the cost and complexity of the higher 
order system increase with increase in order of system. The model order reduction is required for 
minimizing the cost and complexity of the higher order system. The numerator and denominator 
coefficients of the reduced order model are obtained by minimizing integral square error using 


the procedure given in section 3 and are available from the following Tables 2 and 3. 


The four reduced order Kharitonov’s transfer functions are (41-44): 


J 53.4014s +90.1511 
K 79.4364s2 +81.7131s + 89.7877 
3 73,2952 + 90.3032 

T= 


77.256%~ +104.937% + 89.5829 


73 _ 59.9965 + 169.9494 
K 50s? +80s +169.765 


74849994 +169.2678 
K 50s +1208 + 169.9993 


(41) 


(42) 


(43) 


(44) 


Therefore the step responses of the original and reduced order Kharitonov’s transfer functions are 


shown in Figure 2, 3, 4, and 5 respectively. 


Table 2. Typical Parameter used by Differential Evolution for Four Kharitonov’s Transfer Functions 


First Kharitonov Second Kharitonov Third Kharitonov Fourth Kharitonov 
Name Value Value Value Value 
Minimum Maximum Minimum Maximum Minimum Maximum Minimum Maximum 
Population size 50 20 50 20 
CR 0.8 0.8 0.8 0.8 
F 0.5 0.5 0.5 0.5 
Numerator parameter] 50 60 70 80 50 60 80 90 
Numerator parameter2 85 95 90 100 160 170 170 180 
Denominator parameter | 70 80 70 80 50 60 50 60 
Denominator parameter2 80 90 100 110 80 90 110 120 
Denominator parameter3 85 95 80 90 80 90 110 120 
Maximum Generations 100 100 100 100 


Table 3. Comparison of ISE for Four Kharitonov Polynomials 


ISE 
Name of polynomial errr 
First Kharitonov 6.677* 10-4 
Second Kharitonov 0.0031 
Third Kharitonov 0.0075 
Fourth Kharitonov 0.0349 


Minimum 
2.4* 10-3 
0.0048 
0.0089 
0.042 
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5.1. Reduced interval model 
Then the closed loop and open loop interval transfer functions model of reduced are obtained from 


the equations (42)-(45) and are represented as(45-46) 


[53.4014,84.9994)s + [90.151 1, 169.9494 
(45) 


Tk = 
[50, 79.4364)s~ + [80,120]s + [89.5829 169.9993] 
[53.4014,84.9994)s + [90.151 1, 169.9494] 
G, (s,c,d) = 5 
[50,79.4364|s~ + [— 4.9994, 66.598 4s + [— 60.3665, 79.8482] (46) 


Here the open loop reduced order model is an unstable and closed loop reduced order model 1s 
stable. In order to improve the stability of this reduced order model, the PI/PID controller is required. 


5.2. Design of PI controller: 
The transfer function of the PI controller is given by 


K; Ne(s) 
=K eS 
Cpy(s) P + : D. (s) 


Then the closed loop transfer function with a PI controller becomes (47) 


[53.401 4K p,84.9994K p J s* + [90.15 LIK p + 53.4014K ; 16994 94K p + 84.9994K , ] s + 
oe [90..4014K ;,169.9494K, ] 
[50,79.4364]s° + [53.4014K p — 4.9994, 84.99 94K p + 66.598] s? + 


[90.151 IK p + 53.4014K; —60.3665,169.9494K p + 84.9994K, +79,8482] 5 + 


[90.151 1K, ,1699494K ; ] (47) 


Srp Hiiop 


Sic Hesji 





c= { irigimal 
Hs- Prepared 





20.8 | 
= 
5 hi 
AEF | 
arp [— {Iripinal 
al 
E i 4 Ñ iii cl 
Dima (apran } o a = Tiras PO PA i k a 
Figure 2. Step Responses for original higher order Figure 3. Step Responses for original higher 
and its reduced order models for first Kharitonov order and its reduced order models for second 
Kharitonov transfer Function 


transfer function 


From the above equation, the characteristic equation of the closed loop interval system with PI 


controller can be taken as (48) 


[50,79.43 64]s° +[53.4014Kp — 4.9994 84.9994Kp + 66.598}s~ + 


[90.151 Kp + 53.4014K, — 60.3665 169.9494Kp + 84.9994K, +79.8482]s+ [90.151 IK,,169.9494K,]=0 


(48) 


From the above equation, the nominal characteristic equation is (49) 
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Lis) = s? +(0.4759+ 1.0693K p )s? +(0.1505+ 2.0095K p + 1.0693K, )s + 2.0095K, =0 (49) 


The optimal characteristic equation that minimizes the Integral Time Absolute Error (ITAE) 
performance for step input is given by [30] is (50) 


3 2 2 3 
Ê is)=s? +1.750,s? +2.1502s +0 =0 (50) 


Step Hésponst 


; i (Priqinial 
1.2 Proposed | 
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Figure 4. Step Responses for original higher order and its reduced order models for third 
Kharitonovtransfer function 
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Figure 5. Step Responses for original higher order and its reduced order models for fourth 
Kharitonovtransfer Function 


By comparing the above two equations (49) and (50), the nominal PI controller parameters as 


0 _ O _ 
BP ag bs, The nominally designed PI controller is shown in Figure 6. This nominally 


designed PI controller is still unable to stabilize the whole closed-loop system under all possible parameter 
variations. It is therefore desired to redesign a robust PI stabilizing controller for this interval plant starting 
from its nominal one. By applying the necessary and sufficient conditions from the Table 1 of the equations 
for the above 3rd order polynomial of the equation (49), the following set of inequality constraints are 
obtained. In order to make this set of constraints into the feasible closed set, a small positive number ‘es’ is 
introduced into the constraints. Hence the Non-Linear Programming (NLP) optimization problem stated as 
0.2 0,2 
to find the controller parameters, KP and Ñ such that the objective function ia ee HEEN 1S 
minimized, subjected to the following constraints Inequality constraints for proposed method: 


Necessary conditions: 
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—90.1511K,; +£ <0 

—90. 1511 K p — 53.4014 K r + 60.3665 +€ <0 

—53.4014 K p + 4.9994 +€ <0 

Sufficient conditions: 

—(—60.3665+ 53.4014K, + 90.151 IK p j +3*(66.598+ 84.9994K p )* 1699494K, +€ <0 


The linear programming problem consists of two decision variables and four constraints. 


The controller parameters KP and ‘tare restricted to small values by choosing the objective function J 
properly. The purpose of using a small positive number ¢ is to formulate a feasible set closed. In this work, 


the DE technique proposed in section 4 is used to minimize the objective function / to determine the values 


of KP and ‘1 they are tabulated in Table 4. From the Table 4, it is observed that KP and £! values are 
increasing as ‘€’ is increased. This represents the sensitivity of the controller with respect to e. By applying 


the proposed algorithm, then the values of controller parameters are obtained as ŚP = 0.7374 and 


Kt= 0.0012 with £ = 0.5. As per rules of the model, the PI controller designed for reduced order is also 
stabilizes original higher order system. Hence the original higher order system is stabilized more by the PI 
controller from this proposed method for the reduced order system. The closed loop step response of 


the system with a PI controller for both proposed method Giza = 0.7374 and */ = 0.0012) and the method 


given in [23] ("e = 0.5 and É! = 0.1) are shown in Figures 7 and 8 for e = 0.5 respectively. The time 


domain specifications of Figure 7 and 8 are shown in Table 5. The step response comparison of four 
extreme plants with a PI controller obtained by the proposed method and the method given in [23] is shown 
in Figure 9. 


Table 4. Variation of “? and £! for Different Values of ¢ for Proposed Method. 


Controller Set E K, K, 
1 0.5 0.7374 0.0012 
2 1 0.9256 0.0045 
3 1.5 1.1752 0.0094 
4 2 1.1961 0.01367 


Table 5. Time Domain Specifications for Higher Order System with PI Controller for Proposed Method and 
the Method in [23] 


Proposed method Existing method in [23] 
Yopeak Peak S Settling %peak Pak Rise Settling 
Name of e Time time . hootMP Time time TER 
KharitonovPolynomial OVErShgOIME tp(sec) tr(sec) He Gee): “<ONCESIOG tp(sec) tr(sec) 

First 14.6937 3.3096 1.3716 5.6396 32.4147 4.008 1.508 10.808 

Second 25.6613 2.5676 0.4562 5.0282 24.3994 2.488 0.880 8.2771 

Third 23.6792 3.1885 1.3010 8.3247 57.6729 3.820 1.378 26.828 

Fourth 30.7057 1.7341 0.5627 4.7373 42.659 2.166 0.741 8.326 


It has been observed from Figures 7 and 8 that the designed PI controller, which uses the proposed 
stability conditions, robustly stabilizes the plant very quickly when compared to the method given in [23]. 
From Table 5, the designed PI controller stabilizes the plant with lesser time domain parameters 
than the existing method. The proposed method involves four sets of equations for NLP to solve. Thus, 
the proposed method requires less computational complexity than the method given in the literature. Thus, 
the developed PI controller from proposed method is not only stabilizes the reduced order model, but also 
stabilizes the original higher order system. 
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Figure 6. Time response to nominal PI controller Figure 7. Time response of proposed interval system 
design for conventional configuration with PI controller with proposed method 
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Figure 8. Time response of interval system with Figure 9. Time response comparison of interval system 
PI controller using the method [23] with the PI controller using proposed method and 
the method in [23] 


5.3. Design of PID controller: 
The transfer function of the PID controller is given by 


No(s) 
D, (s) 


K 
Cp,(s) = Kp +—+Kps= 
S 
Then the closed loop transfer function with a PID controller becomes (51) 
[53.4014K p 849994Kp ] 5° +[90.1511Kp +53.4014K p,169.9494K p +84.9994K p ] s? + 
7 [90.151 1K p + 53.4014K,,169.9494K p +84.9994K ; ] s+[90..4014K 7 ,169.9494K, ] 


[50+53.4014K p 79.43644 84.9994K p Js? +[ 9OISTIKp +53.4014K p —4.9994, 1699494K p +84.9994K p +66.598] s? + 
[90.151 IK p +53.4014K, —60.3665, 1699494K p + 84.9994K r +79.8452]5+[901511K,,1699494K ; ] (51) 


From the above equation, the characteristic equation of the closed loop interval system with PI 
controller can be taken as (52) 


[50+ 53.4014K p 79.4364+ 84.9994K p ]s? +[ 90.15 11K p + 53.4014K p — 4.9994, 169.9494K p +84.9994K p +66.598] s? + 
[90.1511K p +53.4014K -60.3665, 169.9494K p + 84.9994K, +79.8482] s + [90.151 IK ;,169.9494K; ] (52) 


From the above equation, the nominal characteristic equation is (53) 
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A(s )=(64.7182+ 69.2004K )s? +(13005K + 69.2004K + 30.7996)s* +(9.7407+ 13005K + 69.2004K )s + 13005K, =0 
D D P P I 1 =C (53) 


By comparing the above two equations (50) and (53), the nominal PID control parameters as 


Q_ O_ 0 i i EE 
Bests eh Rat A “Dat. The step response to the nominal PID controller is shown in Figure 


10. This nominally designed PID controller is still unable to stabilize the whole closed-loop system under 
all possible parameter variations. It is therefore desired to redesign a robust PID stabilizing controller for 
this interval plant starting from its nominal one. By applying the necessary and sufficient conditions from 
the Table 1 of the equations for the above 3rd order polynomial of the equation (53), the following set of 
inequality constraints are obtained. In order to make this set of constraints into the feasible closed set, 
a small positive number ‘e’ is introduced into the constraints. Hence the Non-Linear Programming (NLP) 


optimization problem stated as to find the controller parameters, Kp. Kt and *D such that the objective 


7 0 \2 0 \2 ne oe ee ; ; 
function 7 =(4p-4’pY +(Kr-K1) +(Kp-KD) ig minimized, subjected to the following constraints 


Inequality constraints for proposed method: 
Necessary conditions: 


~90.1511K]; +£ <0 

—90.1511K p — 53.4014 Kī +60.3665 +€ <0 

—90.1511K p — 53.4014 K p + 4.9994 +€ <0 

—~53.4014K p —50 +€ <0 

Sufficient conditions: 

~(-60.3665+ 53.4014K, + 90.151 IK p f +3*(66.598+ 169.9494K p + 84.9994K p )* 1699494K, +£ <0 


The linear programming problem consists of three decision variables and five constraints. 
The controller parameters Kp. Kt and ŚP are restricted to small values by choosing the objective function J 
properly. The purpose of using a small positive number e is to formulate a feasible set closed. In this work, 
the DE technique proposed in section 4 1s used to minimize the objective function J to determine the values 
of KP: Stand *®D . They are tabulated in Table 6. From the Table 6, it is observed that Kp, Ki ang Kp 
values are increasing as ‘e’ is increased. This represents the sensitivity of the controller with respect to e. 
By applying the proposed algorithm, then the values of controller parameters are obtained as ŚP = 0.7879, 
Kı = 0.0018 and Kp = 0.1716. As per rules of the model, the PI controller designed for reduced order is also 


stabilizes original higher order system. Hence the original higher order system is stabilized more by 
the PID controller from this proposed method for the reduced order system. The closed loop step responses 


of the system with a PI controller for both proposed method (SP = 0.7879, $1 = 0.0018 and Kp = 0.1716) 
and the method given in [23] CoP =-().5, Ki 0.1 and Kp = 0) are shown in Figures 8 and 11 for ¢ = 0.5 
respectively. The time domain specifications of Figure 8 and 11 are shown in Table 7. The step response 


comparison of four extreme plants with the PID controller obtained by the proposed method and 
the method given in [23] is shown in Figure 11. 


Table 6. Variation of £P ,“/ and £P for Different Values of e for the Proposed Method 


Controller Set E€ K, K; Kp 
1 0.5 0.7898 0.0018 0.1716 
2 1 0.9962 0.0034 0.2538 
2 1.5 1.2308 0.0084 0.2780 
4 2 1.2782 0.0092 0.2954 


IJRA, Vol. 7, No. 4, December 2018: 232 — 250 


IJRA ISSN: 2089-4856 O 247 


Table 7. Time Domain Specifications for Higher Order System with Proposed PID Controller for Proposed 
Method and the Method in [23] 


Proposed method Existing method in [23] 

Name of the % peak Peak ot oe Settling % peak Peak Rise Settling 

. f Rise time . i f . 
Kharitonov overshoot Time irse) time overshoot Time time time 
Polynomial MP tp(sec) ts(sec) MP tp(sec) tr(sec) ts(sec) 
First 9.1431 3.0338 1.4152 5.7966 32.4147 4.008 1.508 10.808 
Second 21.0742 1.3867 0.4000 6.2879 24.3994 2.488 0.880 8.2771 
Third 16.4612 3.1852 1.3010 5.6969 57.6729 3.820 1.378 26.828 
Fourth 19.3460 1.6226 0.5858 4.7784 42.659 2.166 0.741 8.326 


It has been observed from Figures 12 and 8 that the designed PID controller, which uses 
the proposed stability conditions, robustly stabilizes the plant very quickly when compared to the method 
given in [23]. From Table 7, the designed PID controller stabilizes the plant with lesser time domain 
parameters than the existing method [23]. The proposed method involves five sets of equations for NLP to 
solve. Thus, the proposed method requires less computational complexity than the method given in 
the literature. Thus, the developed PID controller from proposed method is not only stabilizes the reduced 
order model, but also stabilizes the original higher order system. 
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Figure 10. Time response to nominal PID controller design for proposed configuration 
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Figure 11. The time response of interval system with PID controller using the method [23] 
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Figure 12. Time response of proposed interval system with PID controller 


6. CONCLUSIONS 

This paper proposes the design of the robust PI/PID controller for higher order interval systems 
using the differential evolution (DE) algorithm. In this paper, a stable reduced interval model is obtained 
from a higher order interval system using the DE. This reduced order interval numerator and denominator 
polynomials are determined by using the Kharitonov’s theorem with the minimization of Integral Squared 
Error using the DE. Then for this lower order interval model, a robust PI/PID controller is designed based 
on the stability conditions for determining robust stability of interval systems with the help of DE. Finally, 
these stability conditions are used to derive a set of inequalities in terms of controller parameters and these 
inequalities are solved to obtain a robust controller with the help of a DE algorithm for an unstable interval 
plant. This designed PI/PID controller for the reduced order model stabilizes the higher order system. 
The proposed PI/PID controller procedure is also applied and demonstrated through a typical numerical 
example. It has been observed from the simulation results that the designed PI/PID controller for 
the reduced order model robustly stabilizes the higher order plant with lesser time domain parameters than 
the existing method. The simulation results are evidence for its robustness in stabilizing the interval process 
using a PI/PID controller. Hence we conclude that the designed PI/PID controller obtained from proposed 
method stabilizes the plant superior than the methods available in the literature. Also, the proposed method 
is simple and involves less computational complexity in comparison with the methods available in 
the literature. 
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